Diffusion on an Ising chain with kinks 
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We count the number of histories between the two degenerate minimum energy configurations of 
the Ising model on a chain, as a function of the length n and the number d of kinks that appear 
above the critical temperature. This is equivalent to count permutations of length n avoiding certain 
subsequences depending on d. We give explicit generating functions and compute the asymptotics. 
The setting considered has a role when describing dynamics induced by quantum Hamiltonians with 
' deconfined quasi-particles. 

I. INTRODUCTION 

£S) . Consider a chain with Ising variables <Ji that can take the values ±1 on n sites. The Ising model is described 
by the following action: S = U^2 i=1 (l — OiO-i+x), where U > 0. This model has two degenerate minimum energy 
configurations. Indeed, both the configurations £1,2 = {c; = =pl, for every i} minimize the energy for this system 
with a value S = 0. Every other configuration can be obtained by flipping the cti variables starting from either Si or 
£2: we see that flipping one variable costs an energy AU. We say that a variable is a kink if the adjacency variables 
have a different value. Flipping an adjacent variable <Ji±\ will not change the total energy of the system. However, if 
we flip a second variable i ± k with k > 2, i.e., we create another kink, then the cost in energy will be 8U. The total 
energy of a certain configuration {/Ji : i = 1, n} will be given by S = AU(d + 1), where d is the number of kinks, 
not including the initial one. 

It is obvious that the two minimum energy configurations Si and £2 are only separated by configurations that 
have energy AU. In fact, we can start from, say £1, pick a site i and flip the variable at from —1 to +1. In this way 
we keep staying in a configuration with energy AU. Then, from there we can reach the configuration £2 by flipping 
only adjacent variables. We say that kinks can propagate freely, that is, without paying energy. The kink is made 
of the walls of a domain with all variables of the same value. In one dimension, making the domain larger does not 
change the domain walls. In statistical mechanics, this implies that the Ising model in one dimension has a critical 
temperature T = 0. At every temperature T > 0, some kinks will be created, and once created, a kink can propagate 
freely and connect the two minima (see, e.g., Q). For this reason, at every T > 0, there is only one minimum energy 
configuration; there are two only when T = 0. 

We call history the permutation on n objects that indicate the sequence of the flips connecting Si to £2- So, for 
instance, a history is given by the string i n ) in which we flip first the variable <Jx, then 02 and so on. As 

we have pointed out, this history requires an energy AU. However, every time we create a kink, we pay an energy 
AU. The energy of a history is therefore the number AU(d + 1), where d is number of kinks created during the entire 
history, without including the first one. It is easy to compute the total number of histories: since we must flip n Ising 
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variable, we can do it in n! possible ways, i.e., the number of permutations of length n. In the present paper we count 
the total number of histories given the parameters n and d. We exhibit generating functions for the sequences and 
observe that these numbers are asymptotically 2 n ~ 2d ~ 1 (d + 1)™. Our technique is the one of generating trees (see, 
e.g., @,|). 

This scenario is related to all those models in condensed matter and statistical mechanics, where two minima are 
separated by excitations that can propagate freely (or quasi-freely). Even if the propagation is not free, i.e., if we add 
a "tension" to the domains, there is a difference between propagating an excitation and creating a new kink 0]. 

Recently such a setting has become important again in the subject of topological quantum computation In this 
context, we have quantum Hamiltonians with dcconfined quasi-particles. The quasi-particles are ends of strings. The 
quantum memory is encoded in the expectation value of large, topological observables in a degenerate ground state. 
A crucial problem is whether in this model quantum memory is stable for temperature higher than absolute zero 
As soon as the temperature is greater than zero, some kinks will be created with some density of probability in time, 
and will start to propagate. If the total time to connect two degenerate ground states will scale with the size of the 
system, the quantum memory will be considered robust, otherwise it will be considered fragile. A precise computation 
will require a knowledge about the number of all possible histories for this process. 

We will deal with a purely combinatorial framework. Let us consider a network modeled by a chain with a finite 
number n of nodes. In the network, the first node is adjacent to second one; the second to the first and the third 
one, and so on. To begin with, let us assume that there is a unique node of the network having a special property. 
For example, the node possesses some information, which has to be diffused to the entire network. When the clock 
ticks, the node passes the information to one its neighbors. The process carries on, and at each time step, a single 
node passes the information. It is clear that it takes exactly n time steps for the information to be transferred to the 
entire network. 

In how many different ways this can be done? Summing over all possible initial nodes, the answer to this question 
is 2™~ 1 . This corresponds to the case with a single kink introduced at the beginning of the process. In Figure [1] 
represented all the 2 3_1 — 4 possible different histories, when the number of nodes is n — 3 (see the vertices 123, 213, 
231 and 321). The time steps of the diffusion are represented from left to right. In this case d — because the only 
kink is the initial one. 

Following the above description, let us now consider a variant in which we introduce extra kinks. Again, at each 
time step a single node passes the information, as we have seen before. We assume however that during the process 
there are exactly d nodes acquiring the information not from an adjacent node. Additionally, these special nodes are 
not pairwise adjacent to each other. The total number of ways in which the information can be propagated to the 
entire network is the number of histories. 

The paper is organized as follows: in Section 2 we introduce the notation and formalize the problem; in Section 3 
we state and prove the results. Section 4 is a brief conclusion. 

II. SET-UP 

The graph P n = (V, E), where V = {l,2,...,n} and edges E = {{1, 2}, {2, 3}, {n - l,n}}, is called n-path. 
We will define below a procedure to assign labels to the vertices of n-paths. The procedure will satisfy the specific 
conditions, giving the diffusion dynamics which we have described in Introduction. For each integer < t < n, we 
define a function f t : V — ► }• This gives the ordered set f t = (/t(l), /t(2), ft{n)), specifying the image of 

the function at time t. When t = and t = n, we assume that fo = (— , — , — ) and /„ = (+, +, +), respectively. 
When t = 1, we chose a vertex v £ V and set fi(v) = +. In other words we designate a vertex to be a kink. At each 
time step t > 2, we chose a vertex w such that 

ft(w)=+ only if f t _ 1 ( z ) = + and {w, z} G E. (1) 

The described rule allows us to transform /o into /„ in exactly n steps. To each one of the ordered sets /o, /i, f n , 
we can naturally associate a permutation 7r of the set of n points [n] = {1, 2, n}. A permutation n is interpreted as 
a word of length n over the alphabet [n], in which every symbol of [n] appears exactly once. Specifically, = j if 
fi{j) = + an d fi-i(j) = —■ Moreover, by following Eq. fl} we have ir(i) = j if there is k < i, such that n(k) = j — 1 
or 7r(fc) = j + 1. 
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Let F n be the set of all permutations of [n] satisfying the above conditions. The fc-th permutation Tt^ k ' S F n is 

associated to the ordered set (/{ , /J , fn)- Notice that k is simply an index to distinguish between maps in F n , 
being F n an unordered set. An an example, the table below includes all elements of F4. The left column represents 

(/i fc \ /2 fc \ •••> /J )■ It is superfluous to include /g fc ' and , since /q' £ ' 1 = (—,—,—,—) and = (+,+,+,+) for 
every k. The right column contains Tr^ k \ The first row is then f± = (+,—,—,—), f% = (+,+,—,—), = 
(+,+,+,-) and vrW = 1234. 



(+,-,-,-),(+,+,-,-),(+,+,+,-) 


1234 


(-,-,+,-),(-,+,+,-),(+,+,+,-) 


3214 


(-,+,-,-),(+,+,-,-),(+,+,+,-) 


2134 


(-,-,+,-),(-,+,+,-),(-,+,+,+) 


3241 


(-,+,-,-),(-,+,+,-),(+,+,+,-) 


2314 


(-,-,+,-),(-,-,+,+),(-,+,+,+) 


3421 


(-,+,-, -),(-,+,+,-),(-,+,+,+) 


2341 


(-,--,+),(-,-,+,+),(-,+,+,+) 


4321 



It is simple to see that F n £ S n if n > 3, where S n denotes the set of all permutations of [n]. Let us consider a 
permutation n £ F4. For instance, let n — 1324. The intermediate steps taking (— , — , — , — ) to (+, +, +, +) induced by 
7r are (+, — , — , — ), (+, — , +, — ) and (+, +, +, — ). Notice that the second step is not in agreement with our procedure. 
This is because 7r(2) = 3 and /2(3) = +, but 7r(l) = 1, while, in order to agree with the procedure, we should have 
7r(l) = 2 or 7r(l) = 4. The permutations in the set S n — F n are then characterized by an extra parameter, which 
we denote by d. This parameter specifies the number of vertices from which the procedure can start independently 
We denote by n^ d ' k ^ the fc-th permutation, with exactly d non-adjacent sequences of +'s, inducing the functions 
fx*' j /a k •••) /n-i ■ When d = 0, we have the above setting. In this case, we simply wrote ir^" 1 instead of 7r(° :fc ). 
The set of the permutations of [n] of the form denoted by F^. There is a clear relation between n and d; 

when n is even, d < n/2; when n is odd, d < (n — l)/2. In the next section, we will find an explicit formula for O^F^. 
It is a simple fact that S n = [J d F d , every n. In the table below, we include the permutations in F% = S4 — F4. Since 
#F 4 = 8, then #F% = 16. 



(+. -• -) 


. (+, +, - 




(+■ 


f,-,+) 


1243 


(-. +. -) 


( + , -, + 




(+,+.+ 


-) 


3124 


(+■ -> -) 


■ (+, -, + 




1 (+ 


+,+,-) 


1324 


(-. +. -) 


( + , -,+ 




(+■-,+ 


+ ) 


3142 


(+. -. -) 


■ (+, -, + 




1 (+ 


-,+,+) 


1342 


(-.-.+.-) 


(-,-.+ 


+ ) 


(+■-,+ 


+ ) 


3412 


(+■ -> -, -) 


■ (+, -, - 


+) 


1 (+ 


+,-.+) 


1423 


(-. -, -■ +) 


(+. -, - 


+ ) 


(+,+, - 


+ ) 


4123 


(+. -> -. -) 


■ (+, - 


+) 


, (+ 


-,+,+) 


1432 


(-. -. -.+) 


(+> -. - 


+ ) 


(+, + 


+ ) 


4132 


(-> +> -. -) 


■ (+,+, - 




, (+ 


+,-,+) 


2143 


(-. -. -.+) 


(-.+. - 


+ ) 


(+, +, - 


+ ) 


4213 


(-, +> -. -) 


. (-, +, - 


+) 


, (+ 


+,-,+) 


2413 


(-. -. -.+) 


(-.+. - 


+ ) 


(-, +, + 


+ ) 


4231 


(-, +, -) 


■ (-, +, - 


+) 




+,+,+) 


2431 


(-■ -, -,+) 


(-. -> + 


+ ) 


(+, + 


+ ) 


4312 



III. ENUMERATION 

A generating tree consists of a root label and a set of succession rules. Following [6J], a generating tree is a rooted 
labeled tree with the property that if v± and «2 are any two vertices with the same label and I is any label, then v\ 
and i>2 have exactly the same number of children with the label To specify a generating tree it therefore suffices 
to specify (i) the label of the root, and (ii) a set of succession rules explaining how to derive from the label of a 
parent the labels of all of its children. For example, a portion of the generating tree Q for the set of all unrestricted 
permutations is illustrated in Figure [1] 

The root is indexed by the permutation of length 1 . A vertex indexed by a permutation 7r of length n has n + 1 
children. These are respectively indexed by the n + 1 permutations that can be obtained by inserting the letter n + 1 
in the word 7Ti7T2 • • • 7r„. Clearly, this tree is isomorphic to a simpler tree, in which the root is labeled by the number 
1, and a node labeled by (n) has exactly n+ 1 children, each one labeled by (n + 1): 

• Root: (1); 



• Rule: (s) ~» (s + l) s+1 . 
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FIG. 1: Generating the permutations according to their length 



In order to enumerate the number of permutations in F^, we define the generating tree G, in which we replace each 
permutation tt = it\tt2 ■ ■ ■ 7r„ with 7Ti = 1 by the vector Vec(7r) = (J, k, r), where ttj = n, tt G F™ and 



as it is shown in Figure [2] 



1, 7T ?l 1 < 7T n _ 1 ; 



(2,0,0) 

7T = 12 




(2,0,1) (3,1,0) (4,1,0) (2,0,1) (3,0,1) (4,0,0) 
tt = 1432 tt = 1342 n = 1324 tt = 1423 tt = 1243 tt = 1234 



FIG. 2: Generating the permutations n — tx\tt2 ■ ■ ■ ir n with 7Ti = 1 according to Vec(-7r) 



Proposition 1 The generating tree Q is given by 

• Roots: (2,0,0) and (1,0,1); 

• Rules: 

(j, k, 0) ^ (1, k + 1, 1) • • • (j, fc + 1, 1)0" + 1, fc, 0) • • • (n + 1, fc, 0); 
(j, fc,l) - (l,fe,l)---(i,fc,l)(i+l,fe,0)---(n+ l,k,0). 

Proof. From the definitions, the permutations 12 and 21 are associated to the vectors (2,0,0) and (1,0,1), respec- 
tively. Let tt — \iT2 • • • 7t„ be a permutation in S n associated to the vector (j, k, 0). By inserting the letters n + 1 
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between position i — 1 and position i, we obtain a new permutation 7r = Itt2 • ■ • 7Tj_i(n + 1)^ • ■ • ir n G SVi+l associ- 
ated to the vector (i, fc + 1, 1), for all i = 2,3,..., j, and the vector (i, k, 0), for all i = j + l,j + 2, . . . , n + 1. Let 
7r = 17T2 • • • 7r„ be a permutation in 5* n associated to the vector (j, fc, 1). By inserting the letters n+ 1 between position 
i — 1 and position i, we obtain a new permutation tt = \tt-i ■ ■ ■ 7Ti_i(n + 1)7^ • • • n n G S'n+i associated to the vector 
(i, k, 1), for all i = 1,2, . . . , j, and the vector (i, fe, 0), for alH = j + 1, j + 2, . . . , n + 1. This completes the proof. ■ 
Let Go(i;u, i>) and G\(t]u,v) be the associated generating functions: 



G P (t; = EE E 



>2 j>2 fc>0 7reS„ :Vcc(ir) = (j,fc,7-) 



= E f "EE^ fc ^ jvfc ' 

n>2 j>2 fc>0 

where r = 0,1. Alternatively, this series counts the vertices of the tree Q by their height and their label (the root 
being at height 2). 

Proposition 2 The generating function for the number of permutations tt = tt\ ■ ■ ■ tt u in S n> n > 2, according to 
the statistics j,k, where TTj = n and tt G F^, is given by Go(t;u,v) + Gi(t;u,v) such that the generating functions 
Go(i;if,u) and Gi(t]u,v) satisfy 

G (t;u,v) = t 2 u 2 + (G (t;u,v) — uGo(tu; l,v)) + (Gi(t;u,v) — uG\{tu; l,v)) 

1 — it 1 — u 

and 

till! til 

G x (i; u, v) = t 2 u + - (G (i; 1, v) - G (t; u, v j) + - (G 1 (t; 1, v) - G x (t; u, v)). 

1 — it 1 — it 

Proof. Since g n ,j,k,r is the number of vertices in the generating tree G, as it was shown in Proposition [1] with label 
(j, k, r), we obtain that 

j'-i 

9n+l,j,k,0 = / ^(ffn,i,fc,0 + 9n,i,k,l), 
z=l 



9n+l,j,k,l — / Affn,-t,fc-l,0 + 9n,i,k,l)- 



1=3 

Let g n +i, r (u,v) = J2j>iY2k>o9 n + 1 <J! k > rU ^ v ^ or r ~ ^>1- Rewriting the above recurrence relations in terms of 
9n,r(u,v), we obtain that 

g n+h0 (u, = E E uJvk E(ff",',fc,o + 9n,i,k,l) J , 

j>lk>l \i=l J 

g n+ i A (u,v) = EE uJ-i;fe ^2(9n,i,k-l,0 + 9n,i,k,l) , 

j>lfc>l \l=j I 



which is equivalent to 



<? n+ l, (w, v) = E E 9n,j,kfiV k {u j+1 + ■■■ + U n+1 ) + E E 9n,3,k,lV k {u j+1 + ■•• + 1i" +1 ), 
j>l fc>0 j>l k>0 

g n+hl (u,v) = J2J29nj,kfiV k+1 {u 1 + ■ ■ ■ + u 3 ) 9nj,k,iv k {u 1 + ■ ■■ +u ] ). 

j>l k>0 j>l k>0 
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Multiplying by t n+1 and summing over n > 2 together with using the initial conditions gi,o(u, v ) — u 2 an d <?2, i( u i v ) = 
it, we obtain that 

Go(t ; u, v) = t 2 u 2 + tj2 E E ^.fc,o* n «*(« j ' +1 + • • • + + * E E E ffn A fc,i* n « fc (« i+1 + ■ ■ • + 

n>2j>lfc>0 n>2j>l/c>0 

= * v + rr^ E E E ^,*,o*V(^ - + ^EEE w.i^ V - 

n>2j>lfe>0 ri>2j>lfc>0 

and 

Gi(t; u, v) = t 2 u + t J2 E E ft'AM^H" 1 + ' ' ' + uj ) + 1 E E E .9»J,M^ V + • • • + u j ) 

n>2j>lk>Q n>2j>lk>0 
n>2j>lfc>0 n>2j>lfc>0 

From the dehnitions of the generating functions Goft;it, v) and Gift; it, u), we get that these functions satisfy the 
following functional equations: 

G ft; u, v) — t 2 u 2 + (G ft; u, v) — uG (tu; 1, v)) + - U (Gift; u, u) - uGiftu; 1, u)), 

1 — U 1 — 14 

and 

twv tu 

Gift; u, ») = t 2 u + (G ft; 1, v) - G ft; u, «)) + (Gift; 1, v) - G x ft; u, »)), 

1—u 1 —u 

as claimed. ■ 

Theorem 3 The generating function 

G (M) = EE# F " Vf 

n>2 d>0 

is given by 



G{tv ^ = J2 4t 2 vT^U(l-(l + 2j)vT^7/ - -/CI -rnr< 



, (1 + vT^5)i+«(l - 2jVT=1^(1 - 2(j + l)tVT^) 2 

Proof. Denote the generating functions Goft; u, v)+Gi(t; u, v) and uGoft; u, v)+G\(t; u, v) by Gft; u, v) and H(t; u, v), 
respectively. Then, Proposition [2] gives that 



(l- t -^0^)G(t;u,v) 



+ ^(1 - u + ut(l - «))G(ut; 1, «) - ^ffft; i, „) + «.« a (n a -i-n 3 *(i-«)) = Q 
and 

+ Gftrf; l,v) - n=^s(l -u-ut(l-v))H(t;l,v) + f^{u - 1 - vu + vu 2 + ut{l - v)) = 0. 

This type of functional equations can be solved systematically using the kernel method (see 0] and references therein). 
In this case, if we substitute u — p — , , ; \. and u = q = - — - 1 .. , we obtain that 



s 1 + t t 2 VT^(2 + VT^ + tVT^) 

Hit; l,v) = G( ; H - ; 
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and 



t 



1 - y/T^vt 1 - y/T^vt 

Hence, the generating function G(t,v) = G(t; l,v) satisfies 
f + y/l-v „, t 



;!,«) + 



ty/T=v) 



(1 - vT=Uf) 2 



G{- 



1 + s/T^vt 1 + y/T^vt 

1 - y/T^V „, t 



,v) + 



t 2 y/T^v{2 + y/Y^i + ty/T^v) 



(1 + y/T^tf 



1 - y/T^t 1 - y/T^t 



v) + 



t 2 y/T^(2 - y/T 



ty/T=v) 



(1 - y/T^t) 2 



which is equivalent to 



G(t,v) 



At 2 y/T 



(f - 2ty/T^) 2 {\ + y/T^U) 



(l-ty/l-V-t(l-v)) + 



-G 



t 



Applying the above functional equation an infinite number of times, we can find an explicit formula for the generating 
function G(t, v), as requested. ■ 

Let h n [v) be the polynomial Yld=o # F n v<i - For " = 2,3,. 



(f + y/T~^) 2 (l-2ty/T^) \l-2ty/T~^ 

s, we can find an e 
10, Theorem [3] gives 



, v 



h 2 (v) 


= 2, 






= 4 + 


2v, 


hi(v) 


= 8 + 


16v, 


h 5 (v) 


= 16 H 


-88w + 16w 2 , 


h 6 (v) 


= 32 H 


-416u + 272v 2 , 


h 7 (v) 


= 64 H 


- 1824v + 2880i; 2 + 272v 3 , 


h s (v) 


= 128 


+ 7680u + 24576w 2 + 7936w 3 , 


h 9 (v) 


= 256 


+ 31616w + 185856u 2 + 137216v 3 + 7936v 4 , 


h 10 (v) 


= 512 


+ 128512u + 1304832i; 2 + 1841152u 3 + 353792u 4 



On the other hand, Theorem [3] can be used to obtain explicit formula for the generating functions h d (t) = 
£ n>2 #F^". For example, Theorem with v = 0, gives h°(t) = ^ = En>2 2n_1 * n : which implies that 
= 2™~ 1 , for all n > 2. More generally, in order to obtain an explicit formula for the generating function 
h d (t), on the basis of the theorem, it is enough to take the derivative of the sum 



2t 2 V j y/T^V 



2j+l 



_ £ (1 - 2jWT=v)(l - 2(j + l)WT=v)(l + y/T^) 

exactly d times and then substitute v = 0. Applying Theorem [3] for d = 0, that is looking for the free coefficient of v 
in the generating function G(t, v) (in this case also j = 0), we get that 

2t 2 

w 1 - 2t 

By considering the coefficient of u 3 , j = 1,2,3, in the generating function G(t,v) (see Theorem [3]) , we obtain an 
explicit formula for the generating function hP{t), and we can state the following result. 
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Corollary 4 We have 



*°(*) = 7^w. ^(*)= 



1-24' w (l-24) 2 (l-44)' 

2 164 5 (l-34) 3 _ 164 7 (17 - 1844 + 6364 2 - 7204 3 ) 



(1 - 24) 3 (1 - 44) 2 (1 - 64)' w (1 - 24) 4 (1 - 44) 3 (1 - 64)2(1- 84) ■ 

This gives an exact formula for #F£, when d = 0, 1, 2, 3 and ra > 1, 3, 5, 7, respectively 

#F° = 2"- 1 ; 

#F„ 1 = 2"- 2 (2"- 1 -n); 

= ^6" - ^4" (n - 1) + ^2" (2n 2 - 4n - l) 

#F " = i4 8 " " el 6 " {n 2) + el 4 " (n2 " 4?1 + 2) ~ ik 2 " (l2n ' " 13n + 2r? + 6) • 

We now focus on the asymptotic behavior of hd- Two sequences a n and b n are said to be asymptotically equivalent 
as n — > oo if linin^oo ^ = 1. We denote this fact by et„ ~ 6„. 

Theorem 5 for any fixed d > 0, 

Proof. Theorem [3] shows that the smallest positive pole of the generating function h d (t) is 4* = 2 {d+i) °^ or der 1. 
Thus, we have ~ Cd2 n (d + 1)". Concerning the constant Cd, we can write 

c d = hm (1 - t/t*) hm f V 4 A g^( 1 - g + 2 ^^ t ~ V - 

' v^odldv d I ^ (l + V T^) 1 +2J(l-2j4 x /l^) 2 (l-2(j + l)4yr^)2 

= hm fl - *1 hm f V 4j!(^t 2 yT^U(l - (1 + 2j)yT^Ut - 4(1 - «)) 

- tig 1 / */* ) „^ 1 ^ dl(i + VT^)1+2j(1 - 2j4\/T^) 2 (l - 2(j + l)4yT^) 2 



= lim lim 
= lim 



44 2 (1 - 4/4*)\/l — v(l - (1 + 2d)VT"^4 - 4(1 - v)) 



t ->t* v^q (i + VT - ^) 1+2d (l - 2dtVT^v) 2 (l - 2(d + l)4\/l — ^)2 

44 2 (1 -t/t*) 2 
£t? 2 1 + M (l-2d4) 2 (l- t/t*) 2 



4(4 



* ^2 



2!+2 d (l - 2c4/4*) 2 
1 



22d+l ' 

as claimed. ■ 



IV. CONCLUSIONS 

In this paper, we have counted the number #F d of histories that allow to connect a ground state of an Ising chain 
with n variables to the other one, having specified the number d of domain walls (or kinks) created in every history. 
We have given the generating functions and studied the asymptotic behavior. 
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The solution of this problem is important in statistical mechanics, to solve the random walk with kinks, and in 
quantum computing, to investigate the robustness of topological quantum memory at finite temperature. There is a 
number of open problems. A first problem would be to count the number of histories to flip the Ising variables along 
a non contractible loop in a torus. A second open problem would be to introduce the possibility that at a given click 
no variable is flipped or that a variables is flipped back. Eventually, the goal would be is to assign probabilities to the 
three different events: flipping an adjacent variable, creating a kink, no flip, and computing the probability of making 
a loop around the torus in a given number of steps. 
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